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ABSTRACT 

The relation between the ratio of infrared (IR) and ultraviolet (UV) flux densities (the infrared excess: 
IRX) and the slope of the UV spectrum (f3) of galaxies plays a fundamental role in the evaluation 
of the dust attenuation of star forming galaxies especially at high redshifts. Many authors, however, 
pointed out that there is a significant dispersion and/or deviation from the originally proposed IRX- 
/? relation depending on sample selection. We reexamined the IRX-/3 relation by measuring the far- 
and near-UV flux densities of the original sample galaxies with GALEX and AKARI imaging data, and 
constructed a revised formula. We found that the newly obtained IRX values were lower than the original 
relation because of the significant underestimation of the UV flux densities of the galaxies, caused by 
the small aperture of IUE, Further, since the original relation was based on IRAS data which covered a 
wavelength range of A = 42-122 /zm, using the data from AKARI which has wider wavelength coverage 
toward longer wavelengths, we obtained an appropriate IRX-/3 relation with total dust emission (TIR): 
log (Ltir/^fuv) = log [l0°- 4 ( 3 06 + 1 - 58 ' 3 ) - l] + 0.22. This new relation is consistent with most of the 
preceding results for samples selected at optical and UV, though there is a significant scatter around it. 
We also found that even the quiescent class of IR galaxies follows this new relation, though luminous 
and ultraluminous IR galaxies distribute completely differently as well known before. 

Subject headings: dust, extinction — galaxies: evolution — galaxies: starburst — infrared: galaxies — 
ultraviolet: galaxies 



1. INTRODUCTION 

The true star formation rate (SFR) of galaxies at var- 
ious redshifts is one of the most fundamental physical 
quantities to understand the formation and evolution of 
galaxies in the Universe. Especially, the SFR density of 
galaxies in a cosmic volume is a very convenient tool to 
explore the galaxy evolution. The conc ept of the cosmic 
SFR d ensity was originally introduced bv lTinslev fc Danlvl 
(1980), and subsequently it has been measured obser- 
vationally at various redshifts up to z ~ 7-8 (e.g.. 



Lilly et al]|l996tlMadau et al.Hl 996. 1991 IConnollv et "ail 
1997t iBouwens et alj|2004 120071 l2009l among others). 

Since the SFR of galaxies is essentially determined by 
measuring observables related to the ultraviolet (UV) flux 
from OB stars, like hydrogen recombination lines or forbid- 
den lines through ionizing UV photons, or UV continuum 
itself. However, the largest uncertainty to determine the 
SFR from these observables is the dust extinction. Since 
active star formation (SF) is always accompanied by dust 
production because of various dust grain formation pro- 
cesses related to the final stage of stellar evolution (e.g., 
Dwek fc Scaldll980l:lpweldll998tlNozawa et al.l [2003ll2010t 
Takeuchi et al.ll2003ll2005cHDrainell2003D . it is always nec- 
essary to deal with the energy "hidden" by dust. The 
absorbed energy by dust is reradiat ed by dust grains at 
far-infrared (FIR) wavelengths (e.g.. ICalzetti et al] |2000: 



iBuat et all 120021) . Now it is shown that the fraction of 
the hidden star forma tion is even dominant at z ~ 1 
( Tak euchi et al.l l2005af) . The most straightforward way 
to remedy the dust extinction problem would be to sum 
up the SFR calculate d from observed (u ncorrected) UV 
and F IR luminosities (|Buat et al. 2007albl; iTakeuchi et al.l 
2010). However, this is not always possible to obtain both 
restframe UV and FIR flux, especially at high redshifts to 
date. 

In general for Local galaxies, when dust grains in galax- 
ies absorb UV light, shorter wavelength flux is more 
strongly attenuated than longer one, which causes redden- 
ing of galaxy spectra^. Then, the strength of reddening 
and absorption fraction of UV are e xpected to be r elated 
in some way. Evaluating this process. iMeurer et al.l (jl999) 
(hereafter M99) proposed a formula to estimate dust at- 
tenuation from the slope of the ultraviolet (UV) spectrum, 
(3. This is based on the tight relation between the infrared 
(IR) to ultraviolet flux ratio (traditionally referred to as 
the infrared excess: IRX) and f3. Since this method does 
not require any other wavelength data to obtain the value 
of attenuation, like IR emission, it was very much appre- 
ciated by researchers of high-redshift galaxies (z ~ 2-5 or 
more). This is because the dust emission of such high- z 
galaxies is usually much more difficult to observe, and their 
UV spectra are ob served at optical wav e lengt hs thanks to 
the redshifts (e.g.. IBouwens et al.ll2007l I2009D . Thus, the 
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IRX-/3 relation is very popularly used in these studies up 
to now. 

However, recently, subsequent studies have shown 



to Local galaxies ( e .g. |Buat et a 


.1 20051: Seibert et all 


2005t Burearella et al. 2005al 


Buat et alJ 2007allbl: 


Cortcsc et al. 2006: Boissier et al. 


20071 Bouauien et all 


2009: Takeuchi et al. 20101 lOverzier et al. 20111). Es- 



pecially, rather quiescently star-forming galaxies tend 
to distribute significantly below the IRX-/3 relation of 
M99. In addition, another extreme categories of starburst 
galaxies, luminous and ultraluminous IR galaxies (LIRGs 
and ULIRGs, res pectively), occupy a region well above 
the M99 relation (iGoldader et al.ll2002HBuat et alj [2005; 
iBurgarella et al.l l2005bl; iTakeuchi et al.1 I2010D . Also, it 
is reported that high-z star-forming galaxies also tend 
to deviate from th e M99 relation fe.g.. ISiana et al.| [2009: 



iReddv et al1l2010[ ). Since the M99's IRX-/3 relation was 
proposed for UV-luminous starburst galaxies, one may 
think it is not surprising that the M99 sample galaxies 
have different property from IR-selected galaxies. In- 
deed, some attempts have been done to explain physi- 
cally the discrepancy and dispersion of t he IRX-/3 relation 
for various categories of galaxies (e.g. Kong et ail 12004 



iBouauien et~aLll2009t IConrov et alj|201olv ~ However, since 
this relation is often applied to high-z galaxies whose na- 
ture is quite unknown, it is reasonable to examine the 
M99 relation according to their original recipe for its firm 
application in future studies. 

The original IRX-/3 relation of M99 was obtained based 
on IRAS and IUE data. Though the latter was an epoch- 
making facility for the UV astronomical studies, it had one 
observational limitation: the aperture of IUE is 10" x 20" 
ellipse with effective circle-equivalent diameter of 14". As 
M99 discussed in their original work, the sample of star- 
burst galaxies in M99 sometimes have large angular diam- 
eters (see Section [2~Tj) . In contrast, the beam size of IRAS 
is quite large, and most of the part of galaxies is covered 
by the IRAS instrument. This difference of apertures of 
IUE and IRAS would cause a problem in the estimation 
of the IRX-/3 relation since it is based on the ratio of the 
two. 

Today, we have a much larger dataset open to public at 
UV bands provided by GALEX. GALEX (Galaxy Evolu- 
tion Explorer) i s a UV astronomical s atellite launched in 
2003 by NASA ()Morrissev et al.ll200 7fl. GALEX has far- 
UV (FUV) and near-UV (NUV) bands, and is performing 
all-sky survey as well as deep surveys in some sky areas. 
In this work, we reexamine the M99 relation by using the 
same galaxy sample as M99, with paying a particular at- 
tention to the aperture problem. For this purpose, we 
use GALEX NUV and FUV images for the M99's orig- 
inal sample galaxies. Another important aspect of using 
GALEX filters for estimating /? of the M99's sample is that 
since the main target of the dust correction by the IRX-/3 
relation is high-z galaxies, for which usually only broad- 
band photometric data are available. Since IUE data are 
spectroscopic, it is crucial to recalibrate photometrically 
estimated UV slope with directly measured one from IUE 



spectroscopic data (see Section 13. 2 .![) . 

It is also interesting and important to examine the IRX- 
j3 relation from the FIR side. For this, we also have 
a large all-sky IR data obtained by AKAR^. AKARI 
is an infrared astronomical satellite launched in 2006 by 
JAXA (Japan Aerospace Exploration Agency). AKARI 
is equipped with two imaging instrumen ts, a near- and 
mid-IR ca mera IRC (lOnaka et al.1 12007D and FIR sur- 
veyor FIS dKawada et all I2007D . and a Fourier spectro- 
graph (|Kawada et al.ll2008l) . In this study, we make use of 
the data obtained by FIS. The FIS has four photometric 
bands, N60 (50-80 /mi), WIDE-S (60-110 /mi), WIDE-L 
(110-180 /mi), and N160 (140-180 /mi). Since AKARI 
has longer wavelength bands than that of IRAS, we can 
estimate the total FIR emission from dust directly without 
any complicated con version formula fe.g. JPale etalll2001l: 
iDale k Heloul [20021 ) . Thus, with GALEX and AKARI 
data, we can obtain an unbiased IRX-/3 relation for the 
M99 sample. 

This paper is organized as follows: in Section [2l we 
describe the data of M99 and measurements newly done 
in this study. Section [31 we first show the result of the 
GALEX-IRAS IRX-/3 relation. Then we examine this 
new relation, and discuss the origin of the discrepancy 
between the M99 relation and recent observations from 
UV side. Then, we examine the IR side with AKARI FIS 
Bright Source Catalog (FISBSC). After considering all the 
measurement-related effects with AKARI diffuse map, we 
propose a new IRX-/3 relation for the UV-selected star- 
bursts with GALEX and AKARI. In Section |H we further 
examine some properties related to the IRX-/3 relation we 
obtained in this work. First we estimate the bolometric 
correction for FIR by IRAS and AKARI flux densities. 
Then, we examined the behavior of the IRX as a function 
of aperture radius but both on the GALEX and AKARI 
images. By this analysis, we can obtain the IRX only from 
the central starburst regions of the original M99 sample. 
After checking these properties, we finally compare our 
new relation to various previous results. Section [5] is de- 
voted to our conclusions. We review the formulation of 
M99 in Appendix[A] In Appendix[Bl we evaluate the aper- 
ture effect in AKARI FISBSC flux densities with AKARI 
FIS diffuse map. 

2. SAMPLE AND DATA ANALYSIS 

2.1. Original data of M99 

First, we descr ibe the original data of M99. 
iKinnev et "all (j 19931) has constructed a UV SED atlas 
of galaxy sample of 143 spiral, irregular, blue compact 
dwarfs, Seyfert 2, and starburst galaxies with IUE spec- 
trum archive. 

M99 selected a subsample of very actively star-forming 
galaxies from this database which consisted of starburst 
nucleus, blue compact, and blue compact dwarfs. Fur- 
ther, in order to reduce the effect of small aperture of 
IUE, M99 have selected galaxies whose surface brightness 
distribution is strongly concentrated in the central small 
region. To construct this sample, they used a criterion 
for an optical diameter £> 2 5 < 4' (£> 2 5 is measured at the 
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25 magarcsec^ 2 isophote). We will revisit this issue in the 
following. 

This further selection resulted in 57 starbursts in the 
IUE sample. M99 used 44 galaxies which were observed 
at IRAS 60 /jm and 100 /mi bands to make the IRX— (3 
plot and obtained the Meurer relation. The sample and 
its property are summarized in Table [TJ 

2.2. New measurement of UV flux by GALEX images 

In order to examine the possible aperture effect of IUE, 
we remeasured the UV flux densities from GALEX images. 
As already mentioned, GALEX has two UV bands, FUV 
(A = 1344-1786 A) and NUV (A = 1771-2831 A). The 
GALEX images are available through MAST^. 

We have made a FUV and NUV photometry of the sam- 
ple galaxies as follows: 

1. Cut out a 30' x 30' square subimage from GALEX 
AIS images around each galaxy 

2. Select a subimage with the largest exposure time 
when multiple observations were available. 

3. Measure FUV and NUV flux densities. The NUV 
image is taken as the reference. 

In this study, it is important to mention that almost all 
of the sources are resolved by GALEX. Then, they very 
often appear as an assembly of small bright patchy re- 
gions, and the GALEX pipeline misidentifies these frag- 
ments as individual objects. This is referred to as shred- 
ding. We must deal with the shredding to obtain sensible 
flux density measurements for nearby extended galaxies. 
For this purpose, we have used an IDL software package 
developed by ourselves. This software performs aperture 
photometry in the NUV sub-image using a set of ellipti- 
cal apertures. Total flux density is calculated within the 
aperture corresponding to the convergence of the growth 
curve. The sky background is measured by combining sev- 
eral individual regions around the source. NUV and FUV 
flux densities are c orrected for Galactic extinction using 
the Schlegel map llSchlegel et all I1998D an d the Galac- 
tic extinction curve of ICardelli et a l. (1989). A detailed 
description of the photome try process can be found in 
llglesias-Paramo et al.l (|2006| ). 

Th is was already u sed fo r previous stud- 

ies (llglesias-Paramo et al.l 120061: iBuat et al.l l2007at 



iTakeuchi et al.T2OT0t ISakurai et al.l 1 2 llh . and its per- 
formance is carefully checked and established. The new 
measurements are tabulated as log Lnuv and log £nuv m 
Tabled] 

3. RESULTS 

3.1. Result with GALEX and IRAS 

We show the result of new measurements for M99's sam- 
ple galaxies to examine the IRX-/3 relation. We quantita- 
tively define the slope of UV spectrum /? as the exponent of 
an approximated spectrum (in flux density per unit wave- 
length interval) 

f\oc\P . (1) 
6 MAST URL: |http://galex.stsci.edu/GR6/| 



The relation between the flux at wavelength A, Fx, and 
the flux density per unit wavelength interval at A, fx is 

Fx = A/a = vfv ■ (2) 
Thus the 1600 A flux Fi6oo is defined as the average flux 
around the central wavelength 1600 A, 

F 1600 = A/ A @1600 A . (3) 
The FIR flux Ffir, is defined by using flux densities /„ at 
IRAS 60 /im and 100 /im as 

F(60) = 2.58 x 10- 40 ^(60 Mm ) , (4) 

F(100) = 1.00 x 10 - 40 /^(i00 /mi) , (5) 
Ffir = 1.26[F(60) + F(100)] , (6) 
where the unit of F(60), F(100), and F FIR is 
erg s" 1 cm~ 2 l, an d the unit of /„ is [erg s _1 cm -2 Hz -1 ] 
Helou et al.lll988T) . In this work, we denote the flux den- 
sity per unit frequency interval at a wavelength A as f v {X). 
We also use flux density in unit wavelength interval, f\. 

The unit of fx is [erg _1 cm _2 s _1 A ] throughout this pa- 
per. In M99, the IR and the UV fluxes were estimated 
from the IRAS and IUE data, respectively. Hence, we first 
use the original IRAS values of M99 for IR flux, and focus 
on the new UV measurements with GALEX. We should 
note that the FIR flux evaluated by eq. (UJ) is only the 
integrated flux over a wavelength range 42-122 /im, and 
not the total IR flux often defined as integrated flux in a 
range 8-1000 /mi. We evaluate the IR flux with AKARI 
later in Section 13.3.21 

With these observables and some empirical assumptions, 
M99 have formulated the relation between IRX1600 and 
Ai 600 as 

log(IRX 1600 ) = log(10°- 4Al600 - 1) + logB , (7) 
BC(1600)„ 

BS BC(FIR) D ;,„ = C °° SL ' < 8 > 
where BC stands for the bolometric correction. M99 
adopted BC(1600)* = 1.66 ± 0.15 for 1600 A and 
BC(FIR) Dus t = 1.4 ± 0.2 for IRAS FIR, which leads 
B = 1.19 ± 0.20. We revisit this issue in Section l+2l 
Thus, we obtain 

log(IRXi 600 ) = log(10°- 4Al60 ° - 1) + 0.076 ± 0.044 . (9) 
The derivation of eqs. (O and © is shown in Appendix IA1 
Here, the value B is for IUE and IRAS, and if we use differ- 
ent data, it will also change. Since we use GALEX FUV, 
we should recalculat e B for GALEX and IRAS. However, 
iSeibert et all (pOOBTl obtained BC(FUV)* ~ 1.68, which 
is very close to that of IUE, yielding practically the same 
value for B. Assuming a linear relation between Aigoo and 
j3 as 

^1600 = ao + aiP , (10) 
namely the general relation between IRX and /3 is ex- 
pressed as 

l og ( = log [io(o-4(«o+«i/0 _ 1 

Performing a least square fit with observed values, M99 
obtained a result 

A woo = 4.43 + 1.99/8. 
Combining eqs. (fT2")l with (jlll) . we obtain 
I Ffir \ 



+ 0.076. (11) 



(12) 



log 



= log 



10 



0.4(4.43+1.99/3) 



0.076. (13) 
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Fig. 1. — The IRX-/3 relation obtained with GALEX and IRAS. Open squares represent the original data of (Mcurcr ct al. 1999: M99), and 
filled circles are the new measurements from GALEX images in this work. A pair of symbols connected by dotted lines represents the same 
galaxy. Typical uncertainty for each measurement is indicated on the lower-left. For the original IUE measurements of M99, the abscissa is 
the UV slope of the IUE spectra of the sample 0, and for the GALEX mea surements, it is the slope defined by the UV color defined by the 
FUV and NUV photometry by GALEX [eq. JTS)], /3 GLX lIKong et al.ll2004D . As for the ordinate, it is the ratio between the FIR flux F F i R 
[eq. II 24 1 11 and the 1600 A flux -F16OO [ ec l- £01 > ''Fill/ I'lCMO (referred to as the IRX) for the original M99 data, and the flux ratio between -FpiR 
and the GALEX FUV-band flux F FUV , F Fm /F Fuv . 



This is the IRX-/3 relation proposed by M99. 

GALEX is a survey-oriented astronomical satellite, and 
its data are mainly images (though GALEX is also 
equipped with a spectrograph). Hence, we need to deter- 
mine the UV spectral slope f3 from FUV and NUV images, 
different from the case of I UE w hich is equipped with a 
spectrograph. iKong et aLI (|2004l ) proposed a method to 
estimate IRX and f3 by using a flux density per unit wave- 
length /fuv and /nuv as follows. 



IRX = log 



FIR. 



/3glx 



-Ffuv, 
log /fuv - log /nuv 



(14) 



(15) 



log Afuv - log Anuv 
Here the effective wavelength of FUV and NUV bands are 
Afuv = 1520 A and Anuv = 2310 A, respectively. 

The IRX-/3 relation newly obtained by these quantities 
is shown in Fig. [1] We also show the original values of M99 
for comparison. Same galaxies are connected with dotted 
lines in Fig. [TJ The new values distribute clearly below 
the original ones of M99 for almost all galaxies. We found 



that since also the new values for j3 are different, the dis- 
tribution is shifted toward lower-right direction. Fitting 
eq. JTTJ to the new GALEX-IRAS IRX-/3 measurements, 
we obtain 

-Ffir 



log 



FUV 



log 



^(0.4(3.45+1.82/3) _ j 



0.076, (16) 



which is different from eq. (fT3"|) (see the dashed line in 
Fig. [T]). The newly obtained relation is very close to 
that of normal sta r-forming galaxie s proposed by previ- 
ous studies (e.g., iBuat et aLl 12005^ | Cortese et al.l [2006; 
Bois sier et alJl2007t iTakeuchi et al. 2010], , among others). 
This implies that the starbursts in M99 sample are hosted 
in a normal galaxies with redder colors. We will revisit 
this issue in Section [4] 

3.2. Examination of the result 

3.2.1. Effect of the filter response function of GALEX 

As we mentioned above, GALEX data are basically im- 
ages, then the photometric data would be affected by the 
filter response function. We should stress that the UV 
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Table 1 

Sample of IMeurer et al"1 (|1999T ) with new measurements. 
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a : fx oc A ^Meurer et al.lH999T ). 

b : ,9glx defined bv lKong etaLI <20olh . 

c : IUE 1600A luminosity IIMeurer et al.|[T999h . 

d : GALEX FUV luminosity. 

c : GALEX NUV luminosity. 

f : IRAS FIR luminosity IIMeurer et al.|[l999l) . 

g : AKARI total luminosity. 

h : Major axis of isophotal ellipse. 

Minor axis of isophotal ellipse. 

Redshift taken from NED. 
k : Morphological type taken from NED. 



spectral slopes /3 were obtained directly from the IUE spec- 
trum, with avoiding the wave length range corresp onding 



trum, with avoiding the wavelength range corresponding 
to the 2170 A bump (see e.g. ICardelli et alJll989j ) of the 
Galactic extinction curve ()Calzetti et al.l ~994f ). In con- 
trast, the GALEXNUV filter covers the bump. This would 
cause a systematic effect on the measurement of IRX and 
p. Since most of the high-z galaxy surveys use photo- 
metric data to evaluate (3 for applying the IRX-/3 relation 
fort dust attenuation, it is crucial to calibrate the related 
quantities carefully. We first examine this effect. 

For this, we convolve the response function with the IUE 
energy spectra of the M99 sample galaxies (see Fig. [2] as 
an example). Let .Rfuv(A) and -Rnuv(A) be the FUV and 
NUV filter response functions, respectively, and denote the 
raw IUE spectra as f\ and the flux density at FUV (NUV) 



as /fuv.AI (/nuv,m) [erg 1 s 1 cm 2 A ], then we have 

f\RFUV(\)d\ 



/fuv 



iii 



FUV 



^?fuv(A)cZA 



FUV 



/a-Rnuv(A)c?A 



/NUV.fil : 



NUV 



Rnuv{X)<I\ 



(17) 



(18) 



NUV 



With eqs. JT7]| and (JT5J), IRX and /3 after transmitting the 
GALEX filters, IRXfu and /3fn, are expressed as 



IRXgi = log 



/ ^FIR 



V-FFUv.ai, 
log /puv,fii - log hvv,m 



(19) 



(20) 



log A FUV - log A NUV 

^FUV.fil = AFUv/FUV,fil • (21) 

We show the comparison of the UV spectral slope /3 and 
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Wavelength [A] 



Fig. 2.— Example of IUE galaxy spectrum in lMeurer et alj 41999ft sample (NGC4861) with GALEX FUV and NUV filter response functions. 
Solid curve represents the observed spectrum of NGC4861 by IUE. Dashed and dotted lines are response functions of GALEX FUV and NUV 
filters, respectively. All are arbitrarily normalized. 
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Fig. 3. — Comparison of the UV spectral slope /3 and the 1600 A flux density before and after transmitting the GALEX filters. Left: 
the ft derived by GALEX FUV and NUV photometry and that directly measured by IUE spectra; Right: the 1600 A flux density through 
the GALEX FUV filter and that directly measured by IUE. Solid lines represent the linear fit to the data, and diagonal dotted lines in both 
panels show y = x. 



the 1600 A flux density GALEX before and after transmit- 
ting the GALEX filter in Fig. [3] For this analysis, we 
restricted to galaxies only with spectra in the full range 
of IUE wavelength coverage (1100-3000 A). In the M99 



sample, there are samples only with spectra at IUE SW 
spectral coverage, and they are not ideal for the purpose 
of this analysis here. As we see in the left panel of Fig. [31 
(3 directly measured by IUE and that defined by GALEX 
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Fig. 4. — Behavior of the IRX-/3 relation with the aperture radius on Gj4Li?.Ximages. The semi-major axis radius of the photometry ellipse 
is set to be 5", 10", 15", 20" from smaller to larger open circles. We also put the converged values of the growth curve of flux given by our 
software by open triangles. 



show a very good agreement. The linear fit gives 

I3 IUE = 1.01/3 GLX ,m - 0.07 (22) 

with the Pearson's correlation coefficient of r = 0.976 and 
the standard deviation is up ~ 0.14. The systematic devi- 
ation is smaller than a p. As for the flux density, as we see 
on the right panel in Fig. [31 the 1600 A flux density mea- 
sured by IUE, Fi 6 oo, tends to be about 10 % larger than 
that transmitted through the GALEX FUV filter, Fpuv,fli- 
The logarithmic linear fit yields 

log fieoo = 0.949 log Ffuv.ai - 0.462 (23) 

with the Pearson's correlation coefficient of r = 0.950 and 
the standard deviation of flux density <jp ~ 0.11 in loga- 
rithmic scale. Again, the systematic deviation is smaller 
than op. Considering the intrinsic scatter in the IRX-/3 
plot in M99, we can safely neglect the systematic differ- 
ences in P and Fpuv or nuv in the following analysis. 

3.2.2. Effect of the aperture size 

In order to examine the effect of aperture photometry 
on GALEX images, we observe the behavior of flux with 
varying the aperture size. In this work, the aperture ra- 
dius 0" of the photometry is defined as an ellipse with 



semi-major axis size 8". This ellipse is determined auto- 
matically by fitting to the isophotal ellipse of a galaxy by 
our photometry software. Then, the flux within a radius of 
10" means the integrated flux in an ellipse with semi-major 
axis size of 10". Figure 2] shows the change of measured 
values on the IRX-/3 plot when we change the GALEX 
aperture radius as 5", 10", 15", 20", and the convergence 
radius (depends on the sample). To see the change more 
clearly, we also show the change of the measured values 
with increasing aperture size with a step of 10" in Fig. [SJ 
In Fig.H we compare GALEX FW and NUV flux den- 
sities and IUE ^1600 flux. Here, the circle-equivalent "ef- 
fective" IUE aperture of 14" is shown as vertical lines. 
This is an approximate of the IUE aperture size which was 
an ellipse of 10" x 20". We set the ellipse shape and posi- 
tion angle to be the same as that used in actual IUE obser- 
vation of the sample galaxies. This IUE effective aperture 
is depicted by vertical lines in both panels of Fig. [6] As 
shown in this figure, FUV fluxes do not start to converge 
around a radius of 10"-20" but continue to increase, and 
finally converge to values approximately 2.5 times larger 
than the 1600 A fluxes. The NUV flux also show a similar 
behavior to that of FUV. We should also note that the 
behavior of FUV and NUV growth curves is not exactly 
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Fig. 5. — Effect of aperture size for the IRX-/3 measurement. These figures are essentially the same as Fig. [4] but we show the change of 
the measured values with increasing aperture size with a step of 10" here. Upper-left: change of the IRX-/3 relation with an aperture radius 
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Fig. 6.— Behavior of GALEX FUV and NUV fluxes normalized by the IUE 1600 A flux. Left: the flux ratio between GALEX FUV and 
IUE 1600 A. Open squares connected by dotted lines represent the flux ratio for each sample galaxy as a function of the aperture radius. 
Solid thick curve is the average of all sample galaxies; Right: same as left panel but for GALEX NUV. A radius 10" is the semi-major axis 
radius of IUE aperture (see main text), which is indicated by vertical dotted lines in both panels. 
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Fig. 7. — Behavior of the difference between /3glx an d Piue as a function of aperture radius of our GALEX photometry. The slope from 
IUE, 0ixjB i s fixed to the original M99 value for each sample galaxy. A radius 10" is the semi-major axis radius of IUE aperture (see main 
text), which is indicated by vertical dashed line. The horizontal solid line means /3glx = Piue- Again the radius 10" is indicated by vertical 
dashed line. 



the same. This causes an aperture dependence of j3. 



To see this effect, we examined the difference of j3 be- 



10 



Takeuchi et al. 



icr 9 r 1 1 1 1 y icr g : 




1 rr" V_ I io~" ki i 

10~" 10~ 10 10~ 9 10~ u 10~ 10 10" 

^fuv t er g s_1 cm -2 ] (Gil de Paz 2007) f NUV [erg s -1 cm -2 ] (Gil de Paz 2007) 



Fig. 8. — Comparison between the measurements of this work and that of Gil de Paz ct al. (2007). Left panel shows the correlation of two 
measurements a t FUV, and right panel is the same but at NUV. In both panels, open squares represent 11 galaxies commonly included both 
in the sample of Gil de Paz ct al. ( 2007) and of this work. Solid lines represent the logarithmic linear fit to the sample. Dotted diagonal lines 
show y = x. 



tween IUE and GALEX. Figure [7] depicts the behavior of 
/?Glx — Piue as a function of aperture radius. We observe 
that the behavior of j3 is not monotonic: some increase 
and others decrease after exceeding 14" aperture, and fi- 
nally start to converge to certain values at much larger 
radius. Within the effective IUE aperture radius, (3glx 
changes significantly for most of the sample and is still 
far from the convergence values. We see a weak tendency 
that the converged values of (3 are larger than that listed 
in M99, i.e., the global UV color of the sample turned out 
to be slightly redder when we measure the integrated flux. 
This can be explained as follows: since the M99 sample 
consists of central starbursts, it may be natural that their 
outer part of the disk has redder colors than the central 
part. Hence, the colors of central part are bluer than the 
global ones. Thus, the values of (3 determined at a radius 
similar to that of IUE are systematically bluer than the 
UV slopes when larger apertures are used. 

We, therefore, conclude that the systematic deviation 
of the IRX-/3 relation from the original one is due to the 
effect of UV flux because of the small IUE aperture. Espe- 
cially, we see there is a galaxy that shows a very prominent 
behavior in Fig. [5] Both at FUV and NUV, this galaxy 
converges to exceptional values much brighter than the 
fluxes at a radius of 10". This galaxy is NGC 2537, which 
has very strong contribution of UV flux from the periph- 
eral regions. Because of this, for NGC 2537, the IUE flux 
is more than an order of magnitude underestimated. 

3.2.3. Verification of the r esult w ith I Gil de Paz et all 

!mM) 

In the previous sections, we have shown that the UV flux 
measurement of M99 was significantly underestimated be- 
cause of the small aperture size of IUE. However, since the 
photometry software we used for t his work has been de- 
veloped originally by ourselves (cf. Ilglesias-Paramo et al.l 



2006), we should examine if there would be any system- 
atic effect caused by the difference in photometry softwares 
w hen we compare our res ults with other ones. 

iGil de Paz et al.l (|2007l ) also used original software to 
perform a photometry of GALEX images. We checked the 
sample of IGil de Paz et all (|2007l ) and found 18 galaxies 
common with M99, among which we hav e 1 1 galaxies with 
meas urements needed for this study (|Gil de Paz et al.l 
120071 see their Table 4). 

Figure HI shows that t he results of this work and that of 
IGil de Pa7et"aT1 (pool agree with each other very well. 
The correlation coefficient and the standard deviation for 
FUV data are r = 0.994 and a F (FVV) = 0.0670 in loga- 
rithmic scale, and the same for NUV data are r = 0.994 
and ctf(NUV) = 0.0630, respectively. Thus, we can 
safely conclude that both measurements are in excellent 
agreement. This means that the aperture effect on 
the original IRX - /3 rela tion was independently proved by 
IGil de Paz etall (l2007lh Now we can safely conclude that 
the reason of the systematic shift of the IRX-/3 relation of 
recent studies and the original one is dominantly caused 
by the small aperture of IUE. Now we have the GALEX- 
IRAS based IRX-/3 relation which is free from the aperture 
effect [eq. ([HI])]. 

3.3. New IRX-P relation with GALEX and AKARI 

Up to now, we concentrated on the examination of UV 
flux measurement for the IRX-/3 relation, and for this pur- 
pose we kept the M99's original values -Ffir from IRAS. 
Here we discuss the IR flux. 

3.3.1. Total IR flux by AKARI 

The M99 relation was constructed to relate the UV slope 
and IRX (i 7 FiR/ / i<i6oo), and they used the FIR luminosity 
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Fig. 9.— New IRX-/3 relation obtained by GALEX and AKARI diffuse map. Measurements with IUE-IRAS, G ALEX-IRAS, and GALEX- 
AKARI are plotted. Diamonds, triangles, and squares represent the IUE-IRAS, GALEX-IRAS, and GALEX- AKARI values. Dot-dashed 
curve represents the original M99 formula [eq. I|13jl ]. dashed curve for GALEX-IRAS fit [eq. I|16| l]. and solid curve for th e new formula wit h 
GALEX- AKARI measurements [eq. pO) ]. Note that for IUE-IRAS a nd GALEX-IRAS, th e IR flux is evaluated by FIR IIHelou et al.llT98gft . 
while for GALEX- AKARI it is evaluated by total IR emission (TIR: ITakeuchi etTaTJE«)T(il) . 



(|Helou et al.lll988l) 

J L F iR = 3.29xl0 _22 x[2.58L^(60^m) + L^(100^m)] [L Q ] . 

(24) 

However, as mentioned above, Lfir contains only the lu- 
minosity in a wavelength range of A = 42 — 122 /im, 
and energy emitted at MIR (5-30 /im) nor submil- 
limeter (300-1000 ^m) is not considered. Hence, the 
FIR luminosity is proved to be a significant underesti- 
mation as a representative of the to tal IR lum i nosity 

~l l2005bri . iDale et al.l (120011) 



Av(WIDE-L) = 0.831 x 10 12 [Hz] , 
logL^ESuu = Ai/( WIDE-S)L V (90 /tin) 



(26) 



-Av(WIDE-L)L u (14D urn) (27) 

(28) 



logLxm = 0.964 log Lf%f RI + 0.814 



(Xtt rI (e.g.. ITakeuchi et 



_jpni_ 

and IDale fc Heloul (|2002l ) proposed formulae to convert 
FIR into TIR luminosity by making use of IRAS flux den- 
sities. 

The .4KA.R/FIS All-Sk y Survey data are pub licly avail- 
able. With AKARI fluxes . ITakeuchi et all (|2010D proposed 
a new simple formula to estimate Ltir precisely as fol- 
lowtQ 



Here Av is the bandwidth of each A KARI FIS fil ter ex- 
pressed in units of frequency (|Hirashita et al.ll2008f) . 

3.3.2. New IRX-/3 relation with GALEX and AKARI 

Combining the GALEX-b&sed UV flux free from the 
IUE aperture effect and the AKARI TIR flux [eq. ([28])] 
covering the whole IR wavelength range, we construct 
a new IRX-/? relation applicable for recent studies with 
AKARI, Spitzer, and Herschel. 

We cross-matched the M99 sample with the database 
of AKARI FIS BSCfl. Among the M99 sample, 48 are de- 
tected by AKARI FIS. Among these 48 galaxies, 40 have 
GALEX FUV and NUV (see Table HJ. With these 40 

7 IHirashita et ahl d200SI) constructed a formula to estimate Ltir by making use of N60, WIDE-S, and WIDE-L flux densities. However, since 
N60 flux density is more noisy than the other wide bands, Takcuchi et al. (2010) found that the formula using two wide bands [WIDE-S and 
WIDE-L: eq. (125 M gives a smaller scatter. Hence, we use this two-band based formula in this work. 

8 DARTS (the AKARI data archive) URL: http://darts.jaxa.jp/astro/akari/cas.html 



Av(WIDE-S) = 1.47 x 10 12 [Hz] 



(25) 
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subsample galaxies of M99, we estimated the new IRX-/3 
relation. 

Again the UV slope of the SED, f3, was obtained by 
eq. ([13), and the IRX was estimated from GALEX FUV 
flux and AKARI TIR flux as follows. 

f -Ftir 



IRX = log 



(29) 



Since the angular resolution of the AKARI diffuse im- 
age map is ~ 30-50" at WIDE-S band, smaller aperture 
radii than 20" are not meaningful. We started an aperture 
radius from 20" which is similar to the IUE aperture ra- 
dius, and then estimated the IRX with 30", 50", and 70" 
radii. We should also take into account the different of the 
PSFs at FIR and FUV. It is pointed out that the synthe- 
sized PSFs of the AKARI diffuse map are more extended 
than the previous estimates, with a significant flux leak 
to the sidelobes (Arimatsu et al., in preparation). Since 
it is still difficult to evaluate this effect precisely at this 
moment, we instead convolved the GALEX FUV images 
with the difference of the angular resolutions of GALEX 
and AKARI, we have convolved the GALEX FUV images 
with the synthesized PSF of AKARI FIS 140 ^m band and 
measured the IRX within the same aperture at FIR and 
FUV bands. The result is shown in Figs. fTUl and [TT1 

It is very clearly seen that the IRX values of central 
regions of the sample are much lower than the relation ob- 
tained from the whole galaxies. This is mainly because of 
the difference of the surface brightness profiles at FUV and 
FIR wavelengths. In general, the FUV light of galaxies is 
much more concentrated than that of FIR. Then, the IRX 
values are well below the convergence values within 15" 
aperture, and then moves up the diagram with increas- 
ing aperture radius, as clearly seen in Fig. 111! Even with 
the 50" aperture radius, the IRX values are still far from 
the convergence values for some galaxies, even if we take 
into account the difference of the PSFs. This implies that 
the starburst regions of the M99's original sample galax- 
ies emit less energy at FIR compared to the whole galaxy 
scale. 

To see this more clearly, we show the fraction of TIR 
flux as a function of aperture radius in Fig. 1121 Though 
the values are largely dispersed, on average the fraction 
contained in 20" aperture radiu£0 is 8.8 ± 3.1 %. Namely, 
though we should keep in mind that the correction factor 
of the flux leak could be large at small aperture radii, at 
least we can claim that less than a few x 10 % of the TIR 
light is associated to the central starburst regions in the 
sample, but comes from outer part of galaxies. 

It will be interesting to simulate the IRX as a function 
of redshift, taking into account a realistic beam size and 
detection limit for current /future observational facilities 
to interpret the high-z IRX-/3 relation. 

4.2. Bolometric correction for the total dust emission 

Since the original IRX-/3 relation was based on IRAS 
observation, we should convert the formula to that of TIR 
and FUV. The assumption of M99 was, as stated above, 
the ratio of TIR/FIR [the bolometric correction of dust 
emission, BC(FIR)D US t: see eq. (J7])] is ~ 1.4. Many au- 
thors used this value to rescale the relation with this fac- 
tor to obtain the IRX defined by log (Ftir/^fuv) (e-g- 
iKong et al.ll2004D . Then, it is important to examine the 
value of BC(FIR)DUSt by new IR data. 

By making use of the AKARI FIS diffuse map, we ex- 
amine the ratio of TIR and IRAS-b&sed FIR. We calcu- 
lated a ratio of TIR/FIR (F Tm /F Fm ), shown in Fig. [H 
In Fig. [T31 the TIR/FIRs among the sample galaxies and 

9 This is practically the smallest meaningful aperture for AKARI diffuse map, as mentioned in the main text. Because of this limited angular 
resolution, we did not try to discuss this issue based on the physical radius of the sample galaxies. 



V -Ffuv , 

The result is presented in Fig. [9] Since the IRAS FIR is 
only a part of the AKARI TIR, the GALEX-IRAS val- 
ues distribute in a lower region than the GALEX- AKARI 
measurements on the IRX-/3 plots. 

However, we should carefully check one issue here. Since 
we tried to be careful about the aperture effect at UV, we 
should also examine if the AKARI BSC flu xes are not 
affecte d by the aperture effect. Importantly, lYuan et ail 
(|201lh examined a sample of Local galaxies extracted from 
that of lTakeuchi et al] (|2010f) , since they are mostly nearby 
resolved galaxies by using AKARI diffuse map. They re- 
measured the IR fluxes with AKARI diffuse images and 
found that the flux of extended galaxies tends to be under- 
estimated in the AKARI FISBSC, which adopted a com- 
plicate d point source extr action procedure. Then, accord- 
ing to lYuan et all (1201 lh . we measured the M99 sample 
with AKARI FIS diffuse map and compared the obtained 
flux with that in the FISBSC (Appendix 0) . Since M99 
tried to remove very extended galaxies in the primary con- 
struction of the sample, the current sample was found to 
be less affected by this effect, but still for some galaxies the 
underestimation was significant. We thus used the fluxes 
from the AKARI diffuse map and used them for the fi- 
nal determination of the IRX-/3 relation. This is shown in 
Fig.! 

For fitting eq. (fTTj) to the data, we should recall that 
this is a formula constructed for IRAS measurement, re- 
flected in the term log B. Since the bolometric correction 
for TIR, BC(TIR) is unity by definition, the last constant 
term in eq. ([TTj) [i.e., eq. (|A10I) ] simplifies to the logarithm 
of the bolometric correction at UV, logBC(FUV) !t ~ 0.22. 
Fitting the eq. (fTT]) to the GALEX- AKARI data with this 
modification yields the following formula. 



log 





= log 


V ^FUV / 





^qO. 4(3. 06+1. 58/3) _ ^ 



0.22. (30) 



This is the new unbiased IRX-/3 relation for UV-selected 
starburst galaxies free from any aperture effect in UV or 
IR, applicable for present-day observational datasets with 
UV images (like GALEX) and total integrated IR flux (like 
that obtained by AKARI). 



4. DISCUSSION 



4.1. 



IRX as a function of aperture radius 

In the above discussion, we have rather focused on the 
aperture effect of the UV photometry on the IRX-/3 re- 
lation. Since we have the AKARI diffuse image map, it 
would be also interesting to explore the radial dependence 
of the IRX values. Here we estimate the IRX-/3 relation 
as a function of aperture radius both on the GALEX and 
AKARI images. For this purpose, we have put the same el- 
liptical aperture homologous to the original IUE one, both 
on these images for each galaxy. 



Reexamination of the IRX-/3 Relation 



13 




Fig. 10. — The IRX-/3 diagram as a function of the aperture size both on GALEX and AKARI images. The same elliptical apertures were 
applied to both of the images. To address the difference of the angular resolutions of GALEX and AKARI, we have convolved the GALEX 
FUV images with the synthesized PSF of AKARI FIS 140 fim band. Since AKARI angular resolution is ~ 50" at longer wavelenghts, 
photometric values with aperture radii smaller than 20" are not shown. Note that the aperture radii are different from that in Fig. [4] they 
are 20", 30", 50", and 70". 



its convergence values are presented. In Fig. [T3l dotted 
lines represent the TIR/FIR ratio, and the symbols repre- 
sent convergence values of TIR/FIR for each sample. The 
thick solid broken line shows the average value. The hori- 
zontal thick line shows the convergence value of TIR/FIR 
at large aperture limit. From this analysis, we found that 
the TIR/FIR ratio converges to 1.44 ±0.574, which is very 
close to th e ratio original l y assumed by M99. By us ing a 
formula of iDale fc HeToul (12002ft . Kong et all (l2004h dis- 
cussed that TIR/FIR is ~ 1.5, which is also quite close to 
our value. 

4.3. Comparison of IRX-/3 relations between various 

samples 

We now compare our new formula with previous stud- 
ies on the IRX-/3 relation with various sa mple selection. 
We su mmarize the comparison in Fig. QH ICortese et all 
( 2006) compiled a sample of optically selected galaxies ob- 
served by GALEX and analyzed the IRX-/3 relation, as 
well as the /ITE-selected UV starbursts. They have con- 
verted the FIR flux to the TIR one by the formula of 

10 URL: http://irsa.ipac.caltech.edu /IRASdocs/hires_over.html. 



IDale et all (|200lD . Since the formula of IDale et ai] (|2001l) 
is defined for a wavelength range of A = 3-1000 /zm, it 
gives a slightly larger value than our formula, but the dif- 
ference is small enough comp ared with the intrinsic scat- 
ter and measurement errors (|Takeuchi et al.ll2Q05bh . To 
describe the average trend, they fit a simple linear func- 
tion instead of the form of eq. pT|) . In their analysis, 
ICortese et al.l (|2006[ ) found that the optically selected sam- 
ple distribute systematically lower than starbursts on the 
IRX-/3 plane. They suggested qualitatively that this is par- 
tially due to the IUE aperture, as well as some physical 
reason related to the star formation activity. Our formula 
passes through the peak of the distribution of their op- 
tical sample, and agrees very well with their linear fit at 
ff >-l. 

iBoissier et al.l (|2007l ) selected 43 nearby late-type spi- 
ral galaxies optically and performed a detailed analysis on 
the dust attenuation. The sample was cross-matched with 
GALEX Atlas of Nearby Galaxies jGil de Paz et aHl2007l) 
for UV and IRAS H I RES image 9^ for IR The lame as 
ICortese et all (|2006l) . IBoissier et all (|2007l ) also used the 
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Fig. 11. — Similar to Fig. [5] but with varying aperture size both on GALEX and AKARI images. The GALEX FUV images are again 
convolved with the synthesized PSF. Since AKARI angular resolution is ~ 50" at longer wavelenghts, we start the aperture radius from 20". 
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Fig. 12. — The fraction of the total IR light within an aperture radius. Dotted lines are the measured fraction for each sample galaxy, and 
the thick broken solid line shows their average. Vertical dashed line represent the semi-major axis of IUE elliptical aperture. The leak of 
flux to the sidelobes of the AKARI PSF is not corrected, but this does not change the conclusion significantly. 



formula of lDale et al.l ()200ll ) to compute the TIR flux from 
IRASHux. For the fitting of the IRX-/3 relation, they used 
a form 



, / -Ftir \ 
l0g liwj 



loe 



iqoo+Oi(FUV-NUV) 



a-2 



(31) 



which has one more free parameter than eq. (|30[) . Their 
sample also distributed below the M99 relation (their 
Fig. 1). Their IRX-/3 relation has a higher IRX at smaller 
/3 and lower IRX at larger j3, but since there are very 
few or no data at f3 < —2, the discrepancy on this side 
is not crucial. At the larger /3 side, actually there is a 
la rge dispersion i n opti cally selected spiral galaxy sample 
in iBoissier et al.1 (120071) . This may be th e reaso n of their 
flat slope. Indeed jMunoz-Mateos et ah (120091) used the 
SINGS galaxy sample (|Kennicutt et al.ll2003l ) for a similar 
analysis and obtained an IRX-/3 relat ion with steeper slope 
at larg e /3 side. The new relation of iMufioz-Mateos et al.l 
(2009) is in a very good agreement with our result at 
(3 > —1.5. Thus, we can safely say that our new IRX-/3 
formula determined by the original UV starburst sample 
of M99, also describes the attenuation of optical sample 
very well. However, we should keep in mind that there 
is a considerable scatter around the relation obtained by 
fitting, as these authors also cautioned. 



Next, we also examine if eq. ([30)1 would apply to the 
IR-selected samples. As we have already seen, the be- 
havior of IR-selected galaxies on the IRX-/3 plane is quite 
complicated. These galaxies have two branches of distri- 
bution, one locates below the M9 9 relation and the other 
spreads well above t he line (e.g., iBuat et ai1l2005l 120101 : 
ITakeuchi et~aT1 1201 Oft. 

iGoldader et al.l "( 20021 ) first discovered that the popula- 
tion distributin g abov e the M99 line is LIRGs or ULIRGs. 
ITakeuchi et al.l (|2010f ) examined this aspect intensively 
and found that the population distributing lower than the 
M99 line is quiescent star-forming galaxies, and confirmed 
that the one above the M99 line is LIRGs or ULIRGs, by 
evaluating the 60-/im luminosity of AKARI FISBSG sam- 
ple. 

IBuat et ahl (|2005l) presented an impressive IRX-/3 plots 
for both GALEXNUV- and IRAS 60 /^m-selected samples. 
The quiescent branch of IR-selected galaxies distribute in 
almost the same way as the UV-selected sample on this 
plane, i.e., below the M99's IRX-/? relation but in parallel 
with it. More active branch galaxies are spread over the 
M99 relation. We see that the former is well described 
by eq. (l3Cfl) without significant bias. However, since the 
distribution of (U) LIRGs are very different from that of 
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Fig. 13. — The ratio between TIR and FIR as a function of the aperture radius. Since FIR values are measured from IRAS, we fix them 
and measured the TIR/FIR (■Ftir/^fir) ratio as a function of aperture radii on AKARI images. Dotted lines represent the TIR/FIR ratio 
as a function of aperture radius, and the symbols represent convergence values of TIR/FIR for each sample. The thick solid broken line shows 
the average value. The horizontal thick line shows the convergence value of TIR/FIR at large aperture limit. Correspondence between the 
data values at 50" radius and convergence values is indicated also by the dotted lines, with a vertical gap to clarify that there is a large radius 
range omitted between them. The leak of flux to the sidelobes of the AKARI PSF is again not corrected. 



quiescent galaxies, our new relation cannot apply to these 
active galaxies. 

As for high-z galaxies, still it is not very easy to ob- 
tain both UV and IR data, bu t gradually d irect measure- 
ments are reported recently. ISiana et all (|2009D showed 
the IRX-/3 measurements for two such objects, MS1512- 
cB58 (z = 2.724) and the Cosmic Eye (z = 3.074), 
both of which are Lyman-break galaxies (LBGs), i.e., UV- 
luminous starbursts. They found that the y locate sig- 
nifican tly below the original M99 relation. iReddv et al.l 
(|2010l ) also discussed the IRX-/3 relation for both UV- and 
IR-luminous galaxies at z ~ 2. They found that their 
"typical" z ~ 2 galaxies do follow the M99 relation (but 
with a large scatter), and ULIRGs at the same redshift 
locate well above it similarly to their Local counterparts. 
Also they argue that young galaxies at these redshifts dis- 
tribute below the M99 line. The young galaxies seem to 
follow the IRX-/? relation expected for the SM C extinction 
curve ([Pettini et al.lll998t ). also suggested by ISiana et al.l 
(2009j). Comparing their results with our new IRX-/3 rela- 
ti on fsee Fig.fT^lj, we would say the z ~ 2 typical galaxies 
of IReddv et al.l (|2010() locate above the new IRX-/3 rela- 



tion, while the young ga l axies follow it. As for lensed z ^ 3 
galaxies in ISiana et al.l (|2009t ). cB58 is approximately on 
the new relation, while the Cosmic Eye may be still be- 
low the relation. It would be very interesting to make 
a statistical analysis when a large sample of galaxies at 
these redshifts will be available, to discuss whether evolu- 
tionary effect really would take place. We should note 
that the IRX-/3 relation of the central starburst regions 
of our sample locates well below the M99's original rela- 
tion and our new result, but more similar to the Pettini's 
relation (see top-left panel of Fig. fTTT) . Since M99 origi- 
nally selected the sample to have a central starbursts, we 
might speculate, for example, that the properties of high- 
z star-forming galaxies could be similar to that of a pure 
starburst region. Further investigation in this direction 
w ill be interes t ing. 

iKong et "ail ([2004) suggested a systematic sequence of 
the IRX-/3 relation from lower to higher along with the 
birthrate parameter b 



b _ SFR(t ) 



(SFR) 



(32) 



where t is a supposed age of a galaxy (jKennicutt et al.l 
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Fig. 14. — Comparison of various IRX-/3 relations in the literature. Thick solid curve: this work; dot-dashed curve: M99's curve converted 
from -Ffir/-Fi600 to Prift/^1600; dashed curve: curve for the M99's sample but with n ew measurements at UV and converted for IR flux, 
i.e., -Ftir/^fuv- long dashed curve: same as -PhrV-Ffuv but with m easurements by Ovcrzicr ct al. (20l j); three- d ot-da shed curve: the 
IRX-/ 3 relation for optica lly selected sample jMufioz- Matcos ct al. 2 00j); solid curve: similar to Munoz-Matcos ct al. (2009) but with older 
data IIBoissier et al.ll2007T) : dotted curve: a (bi sector) linear fit by ICortese et alj i 20061 ): and thin broken solid curve: the IRX-/3 relation from 
the SMC extinction curve l|Pettini et aIJI199Sf ). 



1 1994 ). They claimed that b is related to the attenuation 

as 

A FUV = 3.87 + 1.87Q3qlx + 0-40 log 6) . (33) 
In iKong et all (|2004D . the original M99 relation corre- 
sponds to the extremely bursting group, i.e., b > 3.0, not 
very consistent with the properties of M99's UV-luminous 
starburst sample. We see that our new IRX-/3 relation cor- 
responds to the intermediate b galaxies similar to normally 
star forming galaxies, now consistent with the original M99 

s ample. 

lOverzier et alj (|201l[ ) independently made essentially 
the same analysis with this study but with their own 
GALEX photometry and IRAS -b&sed IR f lux e stima- 
tion with a conversion formula of iDale et all (|2001[ ). with 
Spitzer measurement for a subsample of 12 galaxies Their 
formula is, as we expect, not very different from our result 
[eq. (j3"0)) ]. but slightly deviated upward than ours at large 
/? and high IRX regime. They have measured the UV flux 
up to 2.5-5 times Kron radius of galaxies, their converged 
values are on average twice larger than the original values 
of M99. As we have shown, the IUE flux is underestimated 
with a factor of 2.5 on average, and sometimes the under- 



estimation is of an order of magnitude. He nce, though 
they d iscuss that the quiescent IR galaxies in lBuat et al.l 
(2010]) locate lower than their line, for example, this might 
be due to the UV flux still missed. 

5. CONCLUSION 

The relation between the ratio of infrared (IR) and ul- 
traviolet (UV) flux densities (the infrared excess: IRX) 
and the slope of the UV spectrum (/3) of galaxies plays a 
fundamental role in the evaluation of the dust attenuation 
of star forming galaxies especially at high redshifts. M99 
first introduced a very useful method to estimate the dust 
attenuation by using the relation between the FIR-to-FUV 
flux ratio and the UV slope, the IRX-/? relation. However, 
subsequent studies revealed a dispersion and/or system- 
atic deviation from the original IRX-/? relation of M99, 
depending on sample selection. We reexamined the origi- 
nal IRX-/? relation proposed by M99 by measuring the far- 
and near-UV flux densities of the original sample galaxies 
used by M99 with GALEX imaging data. We summarize 
our conclusion as follows: 

1. The UV flux densities of the galaxies used in M99 
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were significantly underestimated because of the 
small aperture of IUE which was 10" x 20". Con- 
sequently, the IRX values were lower than that in 
M99, which lead the IRX-/? relation shifted down- 
ward on the plot. 

2. The aperture effect affected not only the IRX but 
also /3, because the curve of growth of the flux densi- 
ties is a function of aperture which depends strongly 
on the wavelength. 

3. Correcting the aperture effect, we obtained a new 
IRX-/? relation for UV-selected starburst galaxies 
for the same sample with M99. 

4. The original relation was obtained based on IRAS 
data. Since IRAS bands cover only a wavelength 
range of A = 42-122 /im, the FIR flux calculated 
from IRAS flux does not represent the total IR 
emission from dust. We used data from AKARI 
which has a much wider wavelength coverage espe- 
cially toward longer wavelengths, and obtained an 
appropriate IRX-/3 relation for the data with total 
dust emission. 

5. We also estimated the IRX-/? relation of the sam- 
ple as a function both on the GALEX and AKARI 
images. We found that the IRX-/? relation of the 
central starburst regions of the sample locates sig- 
nificantly below the relation obtained for the whole 
galaxy regions. This means that the extended dust 
emission is important to evaluate various global 
properties of galaxies, though we need to know the 
behavior of the PSFs of AKARI diffuse map more 
precisely to evaluate more quantitatively. 

6. In many previous studies, authors use a con- 
ventional conversion factor (bolometric correction: 
BC) of ~ 1.4 to obtain the total IR flux from dust. 
We examined the BC of dust emission and found 
that indeed the value was 1.44 ± 0.574 on average. 

7. Our new relation is consistent with most of the pre- 
ceding results on the IRX-/3 relation for samples 
selected at optical and UV, though there is a signif- 
icant scatter around it. We also found that even the 
quiescent class of IR galaxies follow this new rela- 
tion, though luminous and ultraluminous IR galax- 
ies distribute completely differently. 



8. The IRX-/? relation for the central starburst regions 
of M99's sample galaxies is very simi lar to the rela- 
tion proposed bv lPettini et al. (1998), based on the 
SMC attenuation curve. This may imply interest- 
ing property of high-z star-forming galaxies. 

In this work, we constructed a proper formula for a cat- 
egory of galaxies which are luminous in UV and actively 
forming stars. We found that this formula also applies 
to various class of moderately star-forming galaxies, like 
optically selected or non-luminous IR selected ones. This 
will be a more firm basis than before when one needs to 
correct dust extinction of high-z galaxies. However, we 
should note the large intrinsic scatter on the IRX-/3 plane, 
especially the (U)LIRG population which does not obey 
the IRX-/? relation at all. Many issues and problems still 
remain to be solved for understanding this relation. 
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APPENDIX 

A. FORMULATION OF IMEURER ET AlI (|1999| ) 

Here, we review the original formulation of the IRX-/3 relation by M99. 

/>oo 

p FLya + / A',0(1 - 10- aAA *')d\' , 

FIR J912A I FIR 



IRX 1600 = = — l£ifi . (Al) 

^1600 F 1600 ,oW- OAAieo ° I W Dust 

where fx'.o is the unattenuated flux density of the emitted spectrum, Ay is the net absorption in magnitudes by dust as 
a function of wavelength, and FLya is the Lya flux. In eq. (|A1|) , it is assumed that the ionizing photons does not affect 
dust heating significantly, namely they are not absorbed directly by dust grains. EE The Lya flux Fhya was produced by 

11 M99 argue that under the starburst population model, typically only < 10 % of the bolometric luminosity of starbursts is emitted below 912 A 
(Leithcrcr & Hcckman 1995), and the fraction of ionizing photons that is directly absorbed by dust grains is estimated to be ~ 25 % llSmith et al.l 
119781 ). However, the direct absorption of Lyman photons by dust grains in ionizing regions has been first discussed byjPctrosian e^aTTi^97^) 
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cascade of photons with wavelengths shorter than 912 A, which finally reach the transition n : 2 — > 1. If the ionized gas 
in Case B, Lya is resonantly scattered repeatedly by hydrogen atoms and finally absorbed by dust. Thus, the numerator 
of eq. (|Aip corresponds to the total flux absorbed by dust, consisting of two terms, the first is the heating by Lya and 
the second is by nonionizing photons. The denominator in eq. (IA1[) represents the observed 1600 A flux density which is 
attenuated by dust. Namely, 

/>OC 

F hoi = F Lya + / . / V)0 (l - 10-° AA >-')dX' , (A2) 



912 A 

^leoo =F 1600,0 10-°- 4Al600 , (A3) 
Pfir \ Flux of dust emission within the wavelength range from 42 to 122 fim (A4) 



F ho \) 



Dust Total flux of dust emission 



Most of the photons absorbed by dust are contribution from massive stars. Then, by assuming that most of the radiation 
from massive stars is UV, we can substitute A\ with ^4i6oo in eq. (|A1I) as 

,0.4A 1600 ^ F ^ + S^2A^''° dX ' ( FfIR 



IRX 160 o = (10 u - 4 ^°°-l)— ; 912A ^ . (A5) 

^1600,0 V -fbol / Dust 

Further, we rearrange eq. (IA5|) as follows. 

IRX 1600 = (10 a4Al6 » - l)B , (A6) 
BC(1600), 

^-BC(FIR) Dust (A7) 
The subscript asterisk means it is related to stellar emission, and subscript Dust means it is for dust emission. Here, 

BC(1600), = ( FLya+ 5 2X/A ' ,0 ' A 1 (A8) 

\ ^1600,0 / 

BC(FIR) Dust = ( . (A9) 

Dust 



/ Fboi \ 
\ -Ffir / 

We assume 



BC(1600)* 

On one hand, BC(1600)„ can be approximated as a c onstant if we fi x the initial mass function (IMF) and a constant 
star formation rate. If we adopt the Salpeter IMF (jSalpeterl Il955| ) with an upper mass limit of 100 M®, we have 
BC(1600)* = 1.66 ± 0.15 (M99). The uncertainty is calculated based on the variation of BC(1600)* depending on the 
burst duration. On the othe r hand, BC(FIR) r ,,,„ ) . is the ratio between the total IR flux and the FIR flux. Empirically 



from the IRAS observation ()Rigopoulou et a l. 1996), M99 adopted a value of BC(FIR)D US t = 1-4 ± 0.2, which leads 



B = 1.19 ± 0.20. Putting them together with eq. (|A6[) . we obtain 

log(IRX 160 o) = log(10°- 4Al6 ™ - 1) + logB 

= log(10 a4j4l6O ° - 1) + 0.076 ± 0.044 . (All) 
which is eq. © in the main text, the basic relation between IRX1600 arL d A\qqq. 

B. EXAMINATION OF THE APERTURE EFFECT IN A KARI FIS BRIGHT SOURCE CATALOG 

The A KARI FISBSC was produced by a dedicated pipeline for point source extraction (jYamamura et al.ll2010l ). How- 
ever, for very extended nearby galaxies, the point source assumption is not valid anymore. To examine this "aperture" 
effect in AKARI FIS flux density, we measured IR flux den sity from AKARI FIS diffuse map. We performed this proce- 
dure by usin g both SExtracto r (|Bertin &: Arnoutil 19961 ) and our own software, and both yielded more or less the same 
result. Since lYuan et al.l (|2011l) used SExtractor, we also present flux density given by SExtractor for the diffuse 
map. We show the comparison between flux densities from FISBSC and diffuse map for M99 sample in Fig. IB11 

Basically these flux densities agree well, because M99 put a constraint to have galaxies with a small angular diameter in 
their sample selection. However, even so, the BSC flux density turned out to be underestimated for non-negligible number 
of galaxies. The horizontal dashed lines represent the average diffe r ence between BSC and diffuse-map flux densities. 
For both WIDE-S and WIDE-L, the average is positive. lYuan et al.l (|2011l) performed a careful comparison of BSC flux 
density with other well-understood database and concluded that this is due to the underestimation of the BSC and not 
the overestimation in diffuse map photometry. Thus, we use the diffuse-map based IR flux densities to have the final 
IRX-/3 relation fitting. 



and extensively by Inoue (2002), who found that the fraction of ionizing photons directly absorbed by dust is significantly larger than these 
values. However, even if all the ionizing photons were absorbed by dust, they would still not drastically affect the conclusions. Hence, the 
formulation of M99 will remain valid. 
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Fig. Bl. — Comparison between the AKARI FISBSC flux density and that measured from the diffuse map at WIDE-S (left) and WIDE-L 
(right). Diagonal lines represent y = x. The lower panels show the logarithmic deviation of data from the diagonal lines. Horizontal dashed 
lines are the average deviation of data at each band. 
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